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Abstract 



We present a study of the inclusive J/ip cross section at y/s = 2.76 TeV and 5.5 TeV. The y/s 
dependence of the cross section, rapidity and transverse momentum distributions arc evaluated 
phenomenologically. Their knowledge is crucial as a reference for the interpretation of A-A and p- 
A J/tp results at the LHC. Our approach is the following: first, we estimate the y/s evolution of the 
Pt-integrated J/tp cross section at mid-rapidity; then, we evaluate the rapidity dependence; finally, 
we study the transverse momentum distribution trend. Whenever possible, both theory driven 
(based on pQCD predictions) and functional form (data driven fits) calculations are discussed. 
Our predictions are compared with the recently obtained results by the ALICE collaboration in 
pp collisions at 2.76 TeV. 
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1. Introduction 

High energy heavy-ion collisions are used to study strongly interacting matter under extreme 
conditions. At sufficiently high collision energies, it is believed that hot and dense deconfined 
matter, commonly referred to as the Quark-Gluon Plasma (QGP), is formed. With the advent 
of a new generation of experiments at the CERN Large Hadron Collider (LHC) [l[ a new energy 
domain is accessible to study the properties of such a state of matter. Many observables have 
been proposed to probe deconfined matter formed in heavy-ion collisions. The production of 
heavy-quarks and quarkonium at LHC energies is promising since their production rate allows 
now careful studies of their characteristics (pt and y distributions, polarization, flow, azimuthal 
correlations) (3,0]. In particular, the production of charm is expected to be large enough to make 
regeneration mechanisms possible. The regeneration of J/ip mesons due to the recombination of 
c and c quarks in later stage of the heavy-ion collision would be a direct probe of the deconfined 
QCD matter jij]. One of the experimental methods to quantify the nuclear medium effects in the 
production of a given observable (Ob) is the measurement of the nuclear modification factor (-R^a) 
in nucleus-nucleus (A-A) collisions, defined as: 

V Ob 

r>Ob _ *AA ( i\ 

~ (N coll ) YO pb W 

where (N co u) is the average number of binary nucleon-nucleon collisions^ and (X^) is the 
invariant yield of the observable Ob in A-A (pp) collisions at a given (same) center-of-mass energy. 
In the absence of nuclear matter effects, the nuclear modification factor should be equal to unity 
for experimental observables commonly called hard probes (large p t particles, jets, heavy- flavour, 
etc). A similar factor B9^, measured in p-A collisions, is crucial in order to disentangle hot and 
cold nuclear matter effects in A-A collisions. 

First results on quarkonium production in proton-proton collisions at y/s = 7 TeV have been 
published by the LHC experiments 0, 0, 0, @] and a successful heavy-ion campaigns took place 
at the LHC in November 2010 and 2011, when Pb-Pb collisions at the energy per nucleon pair 
°f \/ s nn = 2.76 TeV where delivered to the LHC experiments, for an integrated luminosity of 
about 9 /ib -1 . The nominal LHC energies for pp and Pb-Pb collisions are y/s = 14 TeV and 
-y/SNN = 5.5 TeV, respectively. 

In this letter, we present an interpolation of the pt-integrated and differential inclusive J/ip cross 
section in pp collisions to y/s = 2.76 TeV and 5.5 TeV. We establish a procedure to determine the 
J/tjj cross section at arbitrary energy values, which makes possible the calculation of the J/ip nuclear 
modification factors for any colliding system and energy at the LHCd. The interpolation is done in 
three steps. The first step is an energy interpolation of the existing p t - m tegrated J/ip production 
cross section measurements at mid-rapidity (Sec. [3]). The second step is the evaluation of the energy 
evolution of the rapidity distribution (Sec. U]). The calculations are made for the rapidity regions 
where J/ip production is measured at the LHC. The third step is the description of the energy 
evolution of the transverse momentum distribution (Sec. EJ). Our prediction of the J/tp cross 



1 The average number of binary nucleon-nucleon collisions can be estimated by the product of the average nuclear 
overlap function (of the nucleus- nucleus collision) and the inelastic proton-proton cross section [jj. 

2 As an example, the center-of-mass energy per nucleon pair in p-Pb collisions is 3.5-4 TeV at the present LHC 
energy, 7 TeV for pp collisions, and of 8.8 TeV at nominal energy. In the present LHC schedule, a p-Pb run is 
foreseen by the end of 2012. Operation at nominal energy will be commissioned after 2013. 



2 



section and kinematic distributions at y/s = 2.76 TeV are compared with the results recently 
obtained by the ALICE collaboration at mid- and forward-rapidity ET 



2. Available experimental data 

The J/ip interpolation to y/s = 2.76 TeV and 5.5 TeV has been performed by considering: 
PHENIX measurements in proton-proton collisions at y/s = 200 GeV in \y\ < 2.2 CDF 
results in proton-antiproton collisions at y/s = 1.96 TeV in \y\ < 0.6 [l2T ]. and ALICE, CMS 
and LHCb results in proton-proton collisions at y/s = 7 TeV in various rapidity regions covering 
\y\ < 4.5 @, Q) Sj- Table [1] presents a summary of the exploited data. In the following, the 



Experiment 




y range 


Pt range 


(Pt) 




[GeV] 




[GeV] 


[GeV] 


PHENIX [11] 


200 


\y\ < 0.35 


< p t < 9 


1.78 ±0.14 


PHENIX [11] 


200 


1.2 < \y\ < 2.2 


< Pt < 7 


1.62 ±0.10 


CDF [12] 


1960 


\y\ < 0.6 


< p t < 20 


2.51 ±0.09 


ALICE [6J 


7000 


-4.0 < y < -2.5 


< p t < 8 


2.44 ±0.18 


ALICE [6j 


7000 


\y\ < 0.9 


< pt < 7 


2.73 ± 0.44 


CMS [7J 


7000 


1.6 < \y\ < 2.4 


< p t < 16 


2.62 ±0.14 


LHCb [8] 


7000 


2.0 < y < 4.5 


< p t < 11 


2.49 ±0.05 



Table 1: Summary of the measurements considered in this work @, Q, fl EH, E3| • 

statistical and systematic uncertainties of the experimental data have been summed in quadrature. 
We assume J ftp are produced un-polarized. Indeed, if we suppose polarization does not evolve (or 
does slowly) with the interaction energy, its influence on the energy interpolation (cross section 
ratios at different energies or rapidities) should be small for a given experiment. Therefore, we 
neglect the uncertainties related to the polarization. 

2.1. PHENIX results 

production in pp collisions at y/s = 200 GeV has been measured by the PHENIX exper- 
iment at the Relativistic Heavy Ion Collider (RHIC) over a rapidity of —2.2 < y < 2.2 and a 
transverse momentum of pt < 9 GeV/c 111 ]. The inclusive differential cross section as a function 
of rapidity was measured for eleven rapidity bins. At mid-rapidity, we have considered the average 
of the measurements at \y\ < 0.35, that is: 

BR a x da pp /dy\^ <Q35 = 44.30 ± 1.40 (stat.) ± 5.10 (syst.) ± 4.50 (norm.) nb [H]. 

The mean transverse momentum (p t ) has been obtained directly from the experimental data and 

is reported in Table. [TJ 



2.2. CDF results 

The CDF experiment at the Tevatron collider has measured inclusive and prompt J/tjj pro- 
duction in pp collisions at y/s = 1.96 TeV (Run II) over a rapidity of \y\ < 0.6 for all transverse 
momenta from to 20 GeV/c The inclusive J/tp cross section is: 



0-pp(|2/| < 0.6) = 4.08 ± 0.02 (stat.) 



-0.36 
-0.33 



(syst.) lib 



The mean transverse momentum (pt) was obtained from the data and is quoted in Table. [TJ 
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2.3. ALICE results 

The ALICE experiment at the LHC has studied inclusive J ftp production in two rapidity 
regions, \y\ < 0.9 and —4.0 < y < —2.5, in pp collisions at y/s = 7 TeV [jj. The inclusive 
production cross sections are: 
dapp/dy | | y |<o.9 — 5.97 ± 1.02 zib, and 
ddpp/dy |_4.o<j/<-2.5 = 4 -21 ± °- 54 

At forward-rapidity, ALICE [y] has measured the differential cross section in six rapidity bins: 

ddpp/dy |_2.8<y<-2.5 = 5.12 ± 0.99 [A>, 

dapp/dy |-3.i<y<-2.8 = 4.34 ± 0.62 fib, 

do-pp/dy |-3.4<2/<-3.i = 4.64 ± 0.66 /xb, 

ddpp/dy I— 3.7<y<-3.4 = 3.59 ± 0.53 fxb, and 

ddpp/dy |_4.o< y <-3.7 = 3.05 ± 0.54 [A>. 

The (pt) from both the mid- and forward-rapidity measurements are summarized in Table. [TJ 

24. CMS results 

The CMS experiment at the LHC has studied prompt and inclusive J/tp production in \y\ < 2.4 
in pp collisions at y/s = 7 TeV 0]. Only the measurement in the rapidity range 1.6 < \y\ < 2.4 has 
been considered in the present work, since in this rapidity range the CMS acceptance for J/ip goes 
down to p t ~ 0. The pt-integrated inclusive production cross section is: 
BR tl x do-pp/dy |i. 6<2/ <2.4 = 424 ± 50 nb. 

The mean transverse momentum (from data) is quoted in Table. [TJ 

2.5. LHCb results 

The LHCb experiment at the LHC has also measured inclusive and prompt J/tp production in 
the rapidity region of 2.0 < y < 4.5 in pp collisions at y/s = 7 TeV 0]. The inclusive production 
cross section is: 

<7 PP (2.0 < y < 4.5) = 11.65 ± 1.55 ph. 

The mean transverse momentum obtained from data is reported in Table. [TJ The LHCb experiment 

measured the J/ip inclusive cross section for five rapidity bins in the rapidity range 2 < y < 4.5 [jj: 

da P p/dy | 2 .o< y <2.5 = 6.20 ± 0.82 ^ib, (pt) 2 .o<»<2.5 = 2.58 ± 0.19 GeV/c, 

d^p/dy | 2 .5 <3/ <3.o = 5.70 ± 0.69 pb, (p t )2.5< a <3.o = 2.59 ± 0.09 GeV/c, 

da P p/dy | 3 . < y <3.5 = 4.94 ± 0.59 fib, (p t ) 3 .o<v<3.5 = 2.51 ± 0.09 GeV/c, 

d^p/dy | 3 .5< y <4.o = 3.81 ± 0.46 fjb, (pt) 3 .5<j / <4.o = 2.42 ± 0.09 GeV/c, and 

dapp/dy | 4 .o<y<4.5 = 2.64 ± 0.33 fib, {pt) i.o< y <4.5 = 2.32 ± 0.10 GeV/c. 

2.6. Other results not considered in this work 

The J/iJj production in pp collisions was studied also by the STAR experiment at RHIC (y/s = 
200 GeV) Q and ATLAS experiment at LHC (y/s = 7 TeV) 0). However, these measurements 
don't extend down to the low pt region and hence don't provide the pt-integrated cross section 
for their respective rapidity range. Since these cross sections are the basic ingredients of the 
interpolation discussed in this work, these measurements were not considered. Results obtained by 
the DO collaboration ILJI at Tevatron have been discarded here since they are given as a function 
of the J/^ pseudo-raplitjU. SPS results have been excluded since they refer to much lower 
interaction energies. 



3 The J /tp rapidity for a given pseudo-rapidity depends on its transverse momentum: since the other results (and 
our calculation) are given in terms of the rapidity, we can not include this measurement in our work. 

4 



3. Interpolation of 3/tp cross section at mid-rapidity 



The energy evolution of the inclusive J/V> production cross section in proton-proton and proton- 
antiproton collisions at mid-rapidity is shown in Fig. [TJ Here we discuss three approaches to 
determine the pt-integrated mid-rapidity cross section at ^fs = 2.76 and 5.5 TeV. The first method 
is a functional form fitting interpolation, while the other two are based on theoretical model 
predictions, FONLL [lj and LO CEM |3. 



3.1. A functional form fitting approach 

The cross section energy evolution can be described by a functional form, a power-law distri- 
bution: 

f(V~s) = A- (vW)\ (2) 

where A, so and b are free parameters. Although such characterization does not rely on any 
theoretical basis, it reproduces the measured distribution. The results are depicted in FigQ] and 
the interpolation results are reported in Tabled The uncertainties on the cross section interpolation 
to y^s = 2.76 and 7 TeV were estimated by re-doing the fit using the upper and lower limits of 
the PHENIX, CDF and ALICE uncertainties in every possible combination. The envelope of these 
results defines the uncertainty. 





BR U X cr 


2.76 TeV 
5.5 TeV 


2301^ nb 
349^ nb 



Table 2: Mid-rapidity cross section results at 2.76 TeV and 5.5 TeV as interpolated by the power-law function. 



3.2. FONLL based approach 

Another approach consists to use the pQCD energy dependence of charm production cross 
section at mid-rapidity to describe the inclusive J /tp production cross section at mid-rapidity. The 
B feed-down contribution to the inclusive J/ip production cross section is considered via the pQCD 
beauty production cross section, weighted by the decay branching ratio, BR(i? — > J/ip)= BR, as: 



da 
dy 



a x 



do- c 
dy 



+ 2BRx^ 
y =o dy 



y=0 



(3) 



Where a is a free parameter, BR = (1.16 ±0.10) x 10 [16J and the charm and beauty differential 
cross sections are calculated in pQCD at Fixed Order with Next to Leading Log resummation 
The FONLL central predictions were obtained with the CTEQ6.6 [19| parton 



13,11 



FONLL 

distribution functions (PDFs) . The fit of the y/s variation of J /ip inclusive production cross section 
with the function of Eq. [21 using the FONLL central parameters is shown in Fig. [2J 

The FONLL calculation uncertainties were evaluated by varying the PDFs, the quark masses, 
and the renormalization and factorization scales. These uncertainties are defined by the enve- 
lope of the obtained cross section, and are practically described by upper and lower curves. To 
account for the calculation uncertainties in the interpolation procedure, we refit the distribution 
considering these upper and lower curves and recalculate the interpolated cross sections. Finally, 
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0.1 



0.05 



★ Data 

— Power Law fit 



_l_d 



"0 1 2 3 4 5 6 7 

energy (is) [TeV] 

Figure 1: Energy dependence of inclusive 3/tl> production cross section at mid-rapidity as fitted by a power-law 
function: BR U x dcr pp /dy| y=0 ( v /s) = Ax (V s / s o) 6 where A = f 24 ± 9 nb, b = 0.60 ± 0.06 and y/sE = 1 TeV. 



the interpolation is validated by observing how the results deviate when removing, one by one, 
the experimental points with the central FONLL parameterization. The results are reported in 
Table [3l The cross sections evaluated when removing an experimental point to the given energy 



5- 10 



T3 10~ 7 



3, 

"D 10" ! 
X 



DC 

m 



10 s 



• BR x FONLL-b 
□ Data 

A Fit using 1.96 and 7 TeV dat a 
Fit using 0.2 and 7 TeV dat a 
Fit using 0.2 and 1 .96 TeV dat a 
| fit using all data 

I i i i i I i i i i I i i i i I i i 



5 6 7 

energy ( ^s) [TeV] 



Figure 2: Energy dependence of inclusive J/i/j production cross section as fitted by Eq. [3] using charm and beauty 
FONLL calculations with their central parameterization. 



Fit 


Data 




BR U x dcipp 


/dy|j /=0 [nb] per *Js [TeV] 


FONLL 


excluded 


0.2 


1.96 


2.76 


5.5 


7 


central 


none 


39 ± 3 


196 ± 1^ 


241 ± 17 


366 ± 25 


422 ± 28 


lower 


none 


45 ± 3 


201 ± 14 


229 ± 16 


299 ± 20 


332 ± 21 


upper 


none 


41 ±3 


195 ± 1^ 


238 ± 17 


357 ± 24 


409 ± 57 


central 


RHIC 


38 ± 3 


150 ± 15 


235 ± 19 


356 ± 28 


^11 ± 51 


central 


Tevatron 


37 ± 5 


189 ± 22 


232 ± 27 


353 ± 40 


407 ± ^5 


central 


LHC 


41 ±3 


505 ± iff 


253 ± 19 


383 ± 28 


442 ± 32 



Table 3: Inclusive J/^> production cross sections at mid-rapidity obtained with the FONLL based approach. In order 
to distinguish the interpolated values from the other fit values at a given energy, the latter are quoted in italics. 

are in fair agreement (within uncertainties) with the measurements (see Sec. [2]). The extrapolated 
values obtained with the three FONLL curves are compatible, although we note that the fit with 
the lower FONLL uncertainty band tends to give significantly lower results than the others at high 
energies □, this is a source of systematic uncertainty. We remark that the results of the FONLL 
interpolation are in good agreement with those of the functional form fit presented in the previous 
subsection (Sec. 13.1] ). 

3.3. LO CEM based interpolation 

The third approach uses the prompt J/tjj cross section en ergy dependence as computed in the 
Colour Evaporation Model (CEM) at Leading Order (LO) |20j. As in the previous subsection 
(Sec. 13. 2p . the B-decay contribution is estimated via the FONLL beauty predictions. The fit is 
defined such that the normalization of the LO CEM calculation is left as a free parameter, a: 



da 
dy 



dacEM 

a x 



+ 2BRx^ 
y =o dy 



y=0 



y =0 dy 

Different LO PDFs have been considered for this exercis^E Here, we only report the results 



obtained with CTEQ6 [2J and MRST01 [22j. The best fit is obtained with the CTEQ6 PDF set 
and is shown in Fig. [3l Similarly to the other approaches, the interpolation has been validated 
by removing the data points one by one from the fit. The results are summarized in Table |U 
The CTEQ6 and MRST01 derived calculations are consistent with each other. The cross sections 
obtained when removing experimental points are also in good agreement (within uncertainties) with 
the measurements (see Sec. [2]). These results also agree with those of the two previous approaches 
(Sec. O and Sec. 1331). 

3.4- Results 

The values obtained for the inclusive J/ip cross sections at yfs = 2.76 and 5.5 TeV, at mid- 
rapidity, using the functional form, the FONLL and the LO CEM based fits are presented in 
Table [5j together with the x 2 and the number of degrees of freedom of the fit. The fit uncertainties 



4 This effect is understood as related to the discontinuous shape of the FONLL lower-band uncertainties for direct 
charm production [13, [T^j . 

5 It is worth noting that LO CEM based interpolation results using other parton distribution functions are 
compatible within uncertainties with the ones reported here. 
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Figure 3: Energy dependence of inclusive J ftp production cross section as fitted by Eq. 0] with the J/ip LO CEM 
and the beauty FONLL calculations, and the CTEQ6 parton distribution functions. 



Fit 
LO CEM 


Data BR U x da pp /dy\ y=0 [nb] per yfs [TeV] 
excluded 0.2 1.96 2.76 5.5 7 


CTEQ6 
MRSTOl 
CTEQ6 
CTEQ6 
CTEQ6 


none 36 ± 3 198 ± ^ 245 ± 17 370 ± 25 425 ± 23 
none 29 ± 2 ift^ ± ^ 248 ± 18 399 ± 28 470 ± 55 
RHIC 35 ± 3 i9i ± 15 236 ± 19 356 ± 28 410 ± 51 
Tevatron 5£ ± 5 193 ± 23 239 ± 28 362 ± 41 416 ± ^7 
LHC 55 ± 5 205 ± iff 257 ± 19 388 ± 28 446 ± 32 



Table 4: Inclusive J/tp production cross sections at mid-rapidity obtained with the LO CEM based interpolation. 
In order to distinguish the interpolated values from the other fit values at a given energy, the latter are quoted in 
italics. 
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Fit 


a x' 2 (ndf) 


/ \b 

As, ale \ 

FONLL central 
FONLL upper 
FONLL lower 
LO CEM CTEQ6 
LO CEM MRSTOl 


1.5 (1) 

(11.4 ± 0.9) • 10 -3 1.0(2) 
(4.9 ±0.4) • 10" 3 0.57 (2) 
(2.9 ±0.2) • 10" 2 0.06 (2) 
1.94 ±0.15 1.3(2) 
4.7 ±0.4 4.4(2) 


Fit 


BRu x da pp /dy\ y=0 [nb] 
= 2.76 TeV ^ = 5.5 TeV 


A X (y/s/sQ ) 

FONLL central 
FONLL upper 
FONLL lower 
LO CEM CTEQ6 
LO CEM MRSTOl 


230;±P 349±|| 

241 ± 17 366 ± 25 
238 ± 17 357 ± 24 
229 ± 16 299 ± 20 
245 ± 17 370 ± 25 
248 ± 18 399 ± 28 


Result for inclusive J/ip 
Estimate for prompt J/ip 


239i? (model) ± 31 (fit) 350^ (model) ± 51 (fit) 
2151?! (model) ± 28 (fit) 315^ (model) ± 46 (fit) 



Table 5: Extrapolated values of the inclusive J/tp production cross section at mid-rapidity according to the phe- 
nomenological, FONLL and LO CEM fits. The fit and the number of degrees of freedom (ndf) are also reported. 

at y/s = 2.76 (5.5) TeV are of 7% for the LO CEM and the FONLL fits, while they are of 13% 
(15%) for the functional form approach]!. To combine these results we: 

i) discard the LO CEM MRSTOl fit, that has a x 2 /ndf > 2, 

ii) compute the average of the three FONLL values, 

iii) estimate the weighted average of the power-law fit, the FONLL-average (ii) and the LO CEM 
CTEQ6 results. 

We sort the uncertainties of the combined result in two components. One is related to the spread 
of the values obtained with the different approaches, resulting in an uncertainty band covering the 
whole range of values. We will refer to this component of the uncertainty as the model related one. 
The other component is related to the fit uncertainties and includes the influence of the experi- 
mental measurements uncertainties. For this contribution, that we will quote as fit uncertainty, we 
(conservatively) assume the largest of the relative fit uncertainties (the functional form fit one), 
that is of 13% (15%) at 2.76 (5.5) TeV. The combined results for the inclusive J/ip cross section 
at mid-rapidity are quoted in Table 

Finally, an estimate of the prompt J/ip production cross section at mid-rapidity can be reported 
subtracting the B feed-down contribution. The LHCb results at yfs = 7 TeV and 2.0 < y < 4.5 
indicate that the contribution from B-decays is 10% with a 15% relative uncertainty 0]. CDF [l2| 



When dealing with asymmetric uncertainties, we quote as relative error the average between the upper and the 
lower value of the asymmetric uncertainties. 
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measured the B feed-down fraction as a function of pt for p t > 1.25 GeV/c. If we consider that for 
p t < 1.25 GeV/c it is ~ 9%, we obtain a pt-integrated B-decay contribution of 10.5%. Assuming 
that the B-decay contribution to the inclusive J/ip production cross section at 2.76 and 5.5 TeV 
is 10% with a (conservative) uncertainty of 30%, we obtain the prompt J ftp cross section at mid- 
rapidity reported in Table [5j 



4. Rapidity distribution 

The knowledge of the rapidity dependence of J/tp production at 2.76 and 5.5 TeV is crucial 
in order to provide a reference for the measurements performed at large rapidities. Here we focus 
on the estimates at y = 3.25 (ALICE muon spectrometer) and y = 2 (region where the CMS 
acceptance for quarkonia goes down to low pt). We consider and compare the interpolation results 
based on the FONLL and the LO CEM calculations used for the calculation of the mid-rapidity 
cross section (Sec. [3|) with those obtained with a functional form fit approach. For this exercise, 
we consider the measurements at RHIC 11] and at the LHC 0, US]- 



4-1. Interpolation based on theoretical calculations 

The FONLL and LO CEM calculations (see Sec. 13.21 and Sec. I3.3|) provide predictions for the 
production cross section rapidity dependence. In order to test their ability to predict the inclusive 
J/V> rapidity shape at different energies, we fit the measured rapidity distributions at yfs = 200 GeV 
and 7 TeV with a model-derived function. The chosen fit function is similar to the y/s form of 
Eq. [3] and Eq. [J] and is defined as: 



dy 



dc m odel , o d o ^Ofa 
"model X + 2 BR X — 

dy dy 



(5) 



We fix the (model dependent) a m odei parameters to be those obtained from the mid-rapidity fits 
in Secj3] and we let the overall normalization (/3) be the only free parameter of the fit. The fits 
are shown in Figs. H] and [5l and their x 2 values are summarized in Table [6l We observe that 
X 2 /ndf < 2 for all the fits: all the models are compatible with the measured inclusive J/ip rapidity 
distribution at both 200 GeV and 7 TeV. The predictions of these models are expressed here in 
terms of an extrapolation factor from mid- to forward rapidity, defined as: 

/*n,(v) = £(y)/£ ■ (6) 

dy I dy y =o 

The values of f{ orw for the considered models at the energies and rapidities of interest are reported 
in Table [6l The results are a combination of these values. We first compute the average of the 
FONLL estimates and then the average between the FONLL-average and the LO CEM CTEQ6 
value. The uncertainties are defined by the envelope of all the values. The inclusive forward rapidity 
cross sections can then be obtained by multiplying the mid-rapidity cross sections of Table [5] by the 
forward rapidity factors. The uncertainties are split in a correlated (from the mid-rapidity cross 
section) and an uncorrelated (from the /f orw factor) component. The results for the inclusive i/ip 
cross section at forward rapidities are reported on Table El 

Similarly, we can calculate the expected rapidity-differential cross section in a number of ra- 
pidity bins of interest for the LHC experiments. The resulting rapidity distributions at 2.76 and 
5.5 TeV are shown in Figs. 6(a) and |6(b)"| 



7 Fixing the a mo( jci parameter allows us to fix the relative normalization of the prompt and feed-down contributions 
at mid-rapidity. 
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FONLL upper fit 
FONLL lower fit 
LO CEM CTEQ6 fit 

J I I I I I L 



Figure 4: Inclusive J/ip rapidity distribution at ■s/s = 200 GeV [TJ as fitted using the LO CEM and the FONLL 
calculations. The measurement uncertainties are obtained by summing in quadrature the PHENIX systematic and 
statistical uncertainties. 




Figure 5: Inclusive J/ip rapidity distribution at yfs = 7 TeV [1,0, 111 as fitted using the LO CEM and the FONLL 
calculations. The data uncertainties are obtained by summing in quadrature the measurement systematic and 
statistical uncertainties. 
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V~s = 


= 200 GeV 




= 7 TeV 


Fit 


X 2 (ndf) 


P 


X 2 (ndf) 


/? 


FONLL central 


4.1 (10) 


1.10±0.06 


10.9 (11) 


0.87±0.03 


FONLL upper 


5.6 (10) 


0.98±0.05 


13.6 (11) 


0.85±0.03 


FONLL lower 


3.6 (10) 


0.98±0.05 


19.0 (11) 


0.95±0.04 


LO CEM CTEQ6 


3.1 (10) 


1.29±0.07 


15.2 (11) 


0.82±0.03 






/forw(2/) 
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= 2.76 TeV 


v^ = 


5.5 TeV 




y = 2.0 


y = 3.25 


y = 2.0 


y = 3.25 


FONLL central 


0.85 


0.59 


0.88 


0.68 


FONLL upper 
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0.63 


0.90 


0.72 


FONLL lower 


0.91 


0.69 


0.94 


0.80 


LO CEM CTEQ6 


0.87 


0.64 


0.90 


0.73 


Average 


0-87±85S 


64+ - 05 

u - D4 -0.05 


0.901°$ 


71+ - 07 

u -'°-0.05 






BRu x dcr pp 


/<K n 4 [nb] 




2/0 




= 2.76 TeV 




5.5 TeV 


2.0 


208 i 28 (corr.) (uncorr.) 


316 ± 65 (corr.) (uncorr.) 


3.25 


153 ±21 ( 


corr.) _yi (uncorr.) 


256 ± 53 (corr.) til ( uncorr -) 



Table 6: Ratios of the mid- to forward-rapidity 3/if> inclusive cross section scaling as obtained with FONLL and LO 
CEM. The values of the fits x 2 j the number of degrees of freedom (ndf), and the /3 parameters at yfs =200 GeV 
and 7 TeV are reported. The interpolated y-differential production cross sections at yfs = 2.76 and 5.5 TeV are also 
quoted. 
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Figure 6: Inclusive J/ip production rapidity distribution at ^fs =2.76 TeV | (a) | and ^/s = 5.5 TeV |(b)| The ALICE 
measurement at y/s —2.76 TeV are also shown with the total (statistical and systematical) uncertainties. 
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4-2. A functional form fit approach: universal rapidity distribution 

In this subsection we explore a model independent approach to the extrapolation of the rapidity 
distribution, based on the search for an universal energy scaling behaviour in the rapidity distri- 
butions measured at different energies. The advantage of this approach is that there is no a priori 
assumption, and the results only depend on the existing measurements and their uncertainties. To 
compare the y-differential cross sections at 200 GeV and at 7 TeV, these values are renormalized 
by the pt-integrated cross section at y = 0. We draw the obtained values versus: y, y/y max and 
y/ybcammax, where y is the J/tp rapidity, y max = ln^/mj^) and ybeammax = ln(^/s/m p ). The 
resulting distributions are shown in Fig. [71 The similarity between the distributions obtained at 
200 GeV and 7 TeV is remarkable, the most striking case being the scaling with y/y ma x- We then 
choose y/y 

max ** 

scaling variable, and perform global fits to the RHIC and LHC data with suit- 
able functions, the purpose being to define -y/s-independent rapidity distribution functions that can 
be used for the calculations at ^=2.76 and 5.5 TeV. We note that the values used to normalize 
the y-differential cross sections in Fig. [7J their value at y=0, is not an optimal choice since depends 
strongly on a single data point. In the following we will re-evaluate this normalization for each 
functional form. 

The first fit function considered is a Gaussian distribution, we consider that: 

= e -(y/y^?/^ 2 y ( 7 ) 

dy dy y =o 

where o~ y is the width of the Gaussian and da/dy\ y=0 is a normalization factor obtained via 
Gaussian fit separately for the two energies. The result of the fit is shown in Fig.[HJ For comparison, 
the RHIC-only and the LHC-only fit results are also drawn, showing that the fit is stable versus 
excluding data points at a given energy. Also, it is worth noting that uncertainties are reduced 
when considering all data points together. 

We compare the Gaussian fit results to those from a set of polynomial functions, defined as: 

=A-Bx(y/y m ^) n (8) 

dy dy y =o 

where n is an even number (up to 8). The resulting fits are shown in Fig. [9l We deduce from 
the xV n df values, reported in Table that all the functions are compatible with the measured 
y-differential distributions. The measurement uncertainties are not small enough to discriminate 
between quite different rapidity shapes, data driven extrapolations are affected by large uncertain- 
ties. 

The scaling factor f{ OYW (y), defined in Eq.[6l has been evaluated for all the considered functional 
form functions. The results are summarized in Table [71 Though it is observed that polynomial 
functions with n = 6 and n = 8 exhibit a larger x 2 than the ones with n = 2 or n = 4 and 
the Gaussian forms, the experimental measurements do not allow to disregard any fit function. 
We combine these results by computing their weighted average. The uncertainty is given by the 
envelope of their uncertainties. The results are reported in Table [7J For completeness, a similar 
exercise was performed using the beam rapidity scaling (y/?/beam max) instead of the J ftp relative 
rapidity (y/y majx ). The obtained forward rapidity scaling factors differ by less than 2%. Since 
the difference is negligible with respect to the other uncertainties, we will neglect it. Using the 
factors /forw reported here and our estimates of the mid-rapidity cross sections (Table [S]), we 
obtain the inclusive differential cross sections quoted in Table [7J More precise measurements of 
the J/ip rapidity distribution would be needed to disentangle between different fit functions and 
reduce the interpolation uncertainties. 
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proton-proton collisions at 200 GeV (PHENIX [ll|]) and at 7 TeV (ALICE, CMS and LHCb experiments 



15 



Fit 




X 2 /ndf 


A 


B 


Lraussian 






11.6/21 


(*y) 0. 


39 ± 0.02 


Poly-(n= 


2) 




9.2/21 


0.98 ± 0.05 


2.1 ±0.2 


Poly— (n= 


4) 




17.3/21 


0.96 ±0.04 


6.6 ±0.8 


Poly— (n= 


6) 




27.5/21 


0.95 ±0.04 


20 ±3 


Poly— (n= 


8) 




36.3/21 


0.95 ±0.04 


65 ±9 








fforw (y) 




Fit 




v^ = 


2.76 TeV 




5.5 TeV 






y = 2 


y = 3.25 


y = 2 


y = 3.25 


Gaussian 




0.76 ±0.07 


0.48 ±0.07 


0.79 ±0.07 


0.54 ± 0.07 


Poly— (n= 


2) 


0.82 ±0.07 


0.51 ±0.11 


0.85 ±0.07 


0.60 ±0.10 


Poly— (n= 


4) 


0.95 ± 0.05 


0.64 ± 0.09 


0.97 ±0.05 


0.76 ±0.07 


Poly-(n= 


6) 


0.99 ± 0.04 


0.78 ± 0.07 


0.99 ± 0.04 


0.74 ± 0.06 


Poly-(n= 


8) 


1.00 ± 0.04 


0.81 ±0.07 


1.00 ±0.04 


0.91 ±0.05 


Average 




0.951^ 


66+ - 15 


0.96i°;?7 


U -'°-0.21 








BRu x do-pp 
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3.25 




159 ± 22 (corr.) tH (^corr.) 


264 ± 54 (corr.) ±75 (uncorr.) 



Table 7: Ratios of the mid- to forward-rapidity J/ip inclusive cross section scaling as obtained with different functional 
forms. The values of the fits x 2 , the number of degrees of freedom (ndf), and the A and B parameters at <Js —200 GeV 
and 7 TeV are reported. The interpolated y-differential production cross sections at y/s = 2.76 and 5.5 TeV are also 
quoted. 
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Figure 8: Normalised J/?/> production cross section as a function of the relative rapidity y/j/ ma x- Three gaussian fits 
are shown: one was performed by using all data drawn here, other using only RHIC data and the last considering 
only the LHC data. 



4-3. Discussion of the rapidity dependence results 

We observe that the results obtained with the functional form fits (Sec. I4.2p are in agreement 
with those obtained by the interpolation based on the FONLL and LO CEM calculations (Sec. l4.1| ). 
In all cases, the functional form approach involves larger uncertainties, explained by the fact that 
it assumes no a priori knowledge of the rapidity shape. Here, we choose not to combine the results 
obtained with these two approaches. 



5. Transverse momentum distributions 

In this section we study the inclusive J/ip transverse momentum distribution. To compare the 
existing measurements at different energies and rapidity domains, we normalize their p t -integrated 
cross section to unity and plot them together versus the z t variable, defined as: 

Zt=Pt/(pt), (9) 

where the values of {pt) are given in Table [TJ The resulting distributions (see Fig. [TOj) show a 
universal behaviour which can be fitted by the following function: 

1 dV =cx no) 

da/dy dz t dy (l + a 2 z 2 )™' 1 ; 

where c = 2 a 2 (n— 1) follows from the normalization to unity. Requiring (zt) = 1, meaning that (pt) 
computed by the fit function is equal to that of Table[H results in: a = T(3/2) T(n — 3/2) / T(n— 1). 
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(a) n = 2 (b) n = 4 




Figure 9: Inclusive J/')/> production cross section as a function of the J/?/> relative rapidity (y/j/ max ) as fitted by a 
polynomial function, A — B x (y/j/max)" with n = 2 |(a)[ _4|_(b)| 6 |(c)| and 8 |(d)| in proton-proton collisions at 200 GeV 
(PHENIX [nl) and at 7 TeV (ALICE, CMS LHCb [l^ljf 
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Figure 10: Universal z% distribution for all the available experimental data. 



We are then left with only one free parameter n. We do a global fit to all the experimental data 
and obtain: n = 3.9 (see Fig. [T0|) . The systematic uncertainties on the n parameter are evaluated 
by performing three global- like fits excluding one by one the data at 200 GeV, 1.96 TeV and 7 TeV. 
This gives finally: n = 3.9^q'|. To demonstrate the scaling versus Zt variable in more details, we 
plot in Fig. [11] the ratio of the experimental data to the fit function value at each data point. We 
have checked that this fit function describes well also the prompt J/ip data at 7 TeV 0, 0] • Similar 
fit of the Z t sc aling function Eq. (|10p to the T(1S) meson p t spectra measured at 1.8 TeV 23j and 
7 TeV 0, 13] gives n = 3.4. 





Data excluded y/s [TeV] 

0.2 1.96 2.76 5.5 7 


(Pt) 

[GeV/c] 


none 2.31 2.48 
RHIC 2.52 2.51 2.50 
Tevatron 2.17 2.25 2.44 
LHC 2.67 3.01 3.15 


Syst. Unc. 


48% 13% 12% 17% 25% 



Table 8: Inclusive J/ip (p t ) versus y/s and its systematic uncertainty (Syst. Unc). 

The next step to obtain the pt distributions at 2.76 and 5.5 TeV is to estimate the J ftp (pt) 
at these two energies. For this we fit the PHENIX, CDF, CMS, ALICE and LHCb (p t ) data by 
a power-law function of energy. The results are reported in Table [8] To evaluate the calculation 
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systematics, data points at different energies have been removed one by one. The last row of Table[5] 
presents the systematic uncertainty on (p t ), estimated as the difference between the extrapolated 
and measured values for y/s = 200 GeV, 1.96 TeV and 7 TeV, and the spread of the results from 
the global fit and the fits excluding data points for y/s = 2.76 TeV and 5.5 TeV. 

The pt distribution interpolation at 2.76 and 5.5 TeV is done using Eq. (|10p and the correspond- 



ing (pt)- The resulting central distributions and their uncertainty bands are shown in Fig. 12(a) 
and Fig. |12(b)| respectively. 

6. Comparison with the ALICE measurements at ^/s = 2.76 TeV 



The ALICE collaboration did recently measure the inclusive J ftp production at y/s = 2.76 TeV lOJ] 



The inclusive J/tp production cross section at mid-rapidity, BRu x da pp /d?/| y= o = 221±57nb[l0fl, 
is in good agreement with the interpolation value, 239jlf (model) ± 31 (fit) nb, see Table [5j 



The pt-integrated rapidity distribution, see Fig. 6(a) is well described by our interpolation. The 
cross section measurements suggest that the mid- to forward-rapidity factor is: /j^p E (— 3.25 < 
y < -3.0) = 0.6 ± 0.2, and //"° E (-3.5 < y < -3.25) = 0.5 ± 0.2. Comparing these values 
to the interpolation results reported in Tables [6H3 we observe that the theory driven (FONLL 
and LO CEM) and functional form (Gaussian, polynomial-n = 2, polynomial-n = 4) calculations 
reproduce the measured trend. However, the polynomial-n = 6 and polynomial-n = 8 functional 
forms tend to overestimate the measured cross sections at forward-rapidities. The ALICE results 



at y/s = 2.76 TeV 10(] seem then to advocate in favor of neglecting the polynomial-n = 6 and 
polynomial-n = 8 functions of the rapidity interpolation approach for future calculations. If we 
disregard these two functions, the functional form rapidity interpolation factors at y/s = 2.76 TeV 
would be of f forw (y = 3.25) = 0.53±g;^ and f forw (y = 2.0) = 0.87±g;?f , while at y/s = 5.5 TeV 
they would be ff orw {y = 3.25) = 0.64!q];q and ff rw(y = 2.0) = 0.89to?o, to re pl ace those of 



Table [3 The measured transverse momentum distribution, Fig. 12(a) is well reproduced by our 
predictions. 
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Figure 12: Extrapolated normalized p t distribution at 2.76 TeV (a) and 5.5 TeV (b) with the uncertainty bands as 
compared with the ALICE data at s/s = 2.76 TeV. 



7. Conclusions 



A procedure to evaluate the energy dependence of the inclusive (and possibly prompt) 3 ftp in- 
tegrated and differential cross sections in proton proton collisions has been established and used 
to compute the J/i/j cross sections at the same centre of mass energy as the one available (per 
nucleon pair) in Pb-Pb collisions at the LHC (2.76 TeV at present and 5.5 TeV in the near fu- 
ture). The calculation exploits the available experimental measurements in three distinct steps. 
The first step consists in describing the energy dependence of the mid-rapidity p t -integrated cross 
section (Sec. [3]). The expectations from functional form (power-law) and model-driven (FONLL 
and LO CEM) fits have been combined together. The results are summarized in Table The 
interpolation systematics is of 15% (20%) at 2.76 (5.5) TeV. The second step is the determina- 
tion of the scaling factors from mid- to forward rapidities (Sec. 0]). Two approaches have been 
discussed, one based on the theoretically (pQCD) predicted rapidity shape, the other being a func- 
tional form approach, based on the energy scaling of the y/y ma ,x distributions. The results with 
both methods are compatible, although the functional form approach involves larger uncertainties. 
Higher precision measurements are needed in order to better constrain the rapidity distribution. 
The third step is the description of the energy evolution of the transverse momentum distribution 
(Sec. [5]), obtained by assuming an universal scaling of the z t distributions (Eq. I10|) . The predicted 
distributions at i/s = 2.76 and 5.5 TeV are drawn in Figs. 12(a) and 12(b) Our predictions of the 
Pt-integrated and differential cross sections at yfs = 2.76 TeV are in agreement, within interpo- 
lation and experimental uncertainties, with the ALICE measurements at the same energy. These 
comparisons confirm the appropriateness of our interpolation procedure, that can be used to obtain 
a pp reference for the analysis of both A-A and p-A data at the LHC. 



yo 


BR U x da pp /d2/|; n ^ o [nb] 
= 2.76x TeV = 5.5 TeV 




2.0 
3.25 


239+^ ( model ) ± 31 ( fit ) 350+51 (model) ± 51 (fit) 
208 i 28 (corr.)+g (uncorr.) 316 zb 65 (corr.)+y^ (uncorr.) 
153 i 21 (corr.)+i2 (uncorr.) 256 i 53 (corr.)+^g (uncorr.) 



Table 9: Summary of the interpolation results for the inclusive J/ip production cross sections at specific rapidity 
values for \fs = 2.76 and 5.5 TeV. 
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